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The evolution ol hole and spin dynamics in high T c superconductors is studied within the self- 
consistent noncrossing approximation of the t— J model in the small hole density limit. As the doping 
concentration is increased, long-range electron correlations disappear rapidly and the quasiparticle 
energy band becomes considerably narrow. At a small hole density long-range antiferromagnetic 
order is destroyed leading to the inadequacy of spin wave basis approximation near small wave 
vectors. Spin excitations near the antiferromagnetic zone boundary are strongly renormalized and 
damped but they are still well described within spin wave basis approximation. 

PACS numbers: 71.27.+a, 74.25.-q 
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Since the discovery of high temperature superconduc- 
tivity in the copper oxide compounds, spin dynamics in 
those materials has received considerable attention in an 
attempt to understand various anomalous properties of 
the normal state as well as superconductivity itself. In 
the cuprates spin fluctuations are believed as a dom- 
inant interaction channel between charge carriers, be- 
cause strong electron correlations present in the copper 
oxide planes tend to significantly suppress charge fluc- 
tuations. In this respect, recently the strong correlation 
based Hubbard and t — J models have been intensively 
studied by various means p| in order to study spin dy- 
namics in cuprates. At half-filling, the quantum Heisen- 
berg model has exhibited both qualitative and quanti- 
tative agreement with most of the experimental findings 
on undoped (parent) compounds. Spin dynamics away 
from half-filling, however, is still far from understanding, 
although many anomalous properties and even supercon- 
ductivity have been observed in this limit. This lack 
of satisfactory microscopic models away from half-filling 
has lead to introducing various phcnomcnological mod- 
els [HQ]. In the present paper we extend the single hole 
problem studied previously by various authors [|| , to the 
case with a finite but small density of doped holes within 
the well-established self-consistent noncrossing approxi- 
mation. As a result of the self-consistent nature of the 
calculation, hole dynamics is also studied on the equal 
footing. Our aim here is to examine the evolution of hole 
and spin dynamics in the small doping regime from the 
insulating (undoped) side. 

The t — J Hamiltonian reads 
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the same site. After Fourier and Bogoliubov transforma- 
tions to Eq. (jl|), the t — J model Hamiltonian becomes 
within the spin wave basis approximation j^j, 
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where and a$ denote a hole annihilation operator with 
momentum k and a spin wave annihilation operator with 
momentum q, respectively. - and Nj: ~ are the bare 
interaction vertices, and Ox is the spin wave energy with 
momentum q. See Refs. H for details, fj, is the chemi- 
cal potential controlling the doping concentration. The 
hopping term of the t — J model transforms into a cou- 
pling term proportional to t due to the on-site no double 
occupancy constraint. Since spins are ordered and thus 
spin waves are well defined excitations within the spin- 
spin correlation length, linear spin wave approximation 
is expected to be good for a small density limit of doped 
holes which the present paper is focused on. 

In this paper we employ the self-consistent noncross- 
ing approximation in the electron self-energy and spin 
excitation polarization functions. The vertex correc- 
tion neglected in this approximation was shown by Sher- 
man H to be of order of the doping concentration so 
that it is quite a reliable approximation for the present 
study. Within the self-consistent noncrossing approxima- 
tion, the hole self-energy and the 2x2 matrix spin po- 
larization functions are given in terms of a 2 x 2 matrix 
spin excitation propagator D(q) and a two component 
interaction vertex M, 



where J is the spin exchange constant. The cf a and Ci j(T 
are electron creation and annihilation operators acting in 
the reduced Hilbcrt space without double occupancy at 
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where T is the temperature and M is defined by 
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Here fc and q are compact notations for (k,ik n ) and 
(q,iq n ), respectively, where k n and g„ are Fermionic and 
Bosonic Matsubara frequencies. Through the hole self- 
energy and spin polarizations calculated above, we obtain 
the interacting hole and spin excitation Green's functions 
after taking the analytic continuation to get retarded 
functions, 
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where the inverse noninteracting spin wave propagator 
[D 0,R (q, y)] _1 is given as 



[D 0, (q, v)] = 
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Because of its matrix structure, D(q) contains four 
terms associated with two diagonal (normal) and two off- 
diagonal (anomalous) components. Here we take spin 
wave excitations as the zeroth order or unperturbed 
states as seen in Eq. (^) and the present calculation can 
deal with the renormalization and damping of the spin 
excitation modes by taking into account spin polariza- 
tions, U R (q, v), in the spin excitation Green's function. 
Furthermore, a possible inadequacy of the spin wave basis 
approximation in some region of the Brillouin zone may 
be gauged by the appearance of the off-diagonal compo- 
nents of the matrix spin excitation propagator, which will 
be discussed later. It should be noted that in Eqs. (||) 
and (H we did not attach any artificial damping factor 
in the denominator of both the hole and spin excitation 
Green's functions in order to obtain the correct hole spec- 
tral functions without smearing of any fine structures. 
Eqs. (§) and (§) are iterated until required convergency 
is reached. In this way both the hole and spin exci- 
tation Green's functions are computed self-consistently. 
Throughout the paper, J/t — 0.4 and the zero temper- 
ature limit (T = 0) are used and a lattice of 20 x 20 is 
used to discretize momentum space. Recently Sherman 
used similar equations to study the doping dependence 
of hole and spin dynamics in cuprates. 

Figures [l] show the evolution of the hole spectral func- 
tions at (7r/2,7r/2) and at (0,0) for doping concentra- 
tions 0.0, 0.95, 2.93 and 5.54%. Due to the finite size 
of the lattice (20 x 20) in the Brillouin zone we could 
not vary the hole density in a continuous manner but 
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FIG. 1. Hole spectral functions at (7r/2,7r/2) and at (0,0) 
for four different doping concentrations, (a), (b), (c) and 
(d) are the hole spectral functions at (n/2,n/2) for 0.0, 0.95, 
2.93 and 5.54%, respectively, (e), (f), (g) and (h) are the 
hole spectral functions at (0, 0) for 0.0, 0.95, 2.93 and 5.54%, 
respectively, where the vertical axis is multiplied by three to 
that for (a')-(d). 



instead in a discrete way as indicated above. The most 
striking change of the spectral functions at (ir/2,ir/2) 
by doping is the significant suppression of string excita- 
tions. For 5.54% hole doping (Fig. 0(d)) these multiple 
magnon excitations become strongly damped and are ab- 
sorbed into a featureless continuous background. Since 
higher string excitations correspond to longer-range cor- 
relations, they are more susceptible to doping, as seen 
in the figure. The strong suppression of the string exci- 
tations for 5.54% doping concentration implies that the 
effective hopping rate i e // already becomes comparable 
to J, leading to the considerable loss of strong correla- 
tions even at this low doping. 

The hole spectral functions at (0, 0) are shown in 
Figs. Jj](e)-(h) for the same doping concentrations as in 
Figs.|[(a)-(d). Although the spectral weight of the lowest 
peak at (0, 0) is relatively smaller than that at (vr/2, tt/2) 
for x = 0.0% a well-defined quasiparticle state with an 
infinite lifetime exists near the bottom of the spectrum 
(Fig. |l|(e)) at (0,0) point. For only a few percent of 
doping concentration, however, the quasiparticle state 
at this momentum is completely destroyed and the in- 
coherent part becomes more dominating. This can be 
easily understood by examining the real and imaginary 
parts of the inverse hole Green's function as exhibited 
in Figs. |^. While the pole condition (vanishing real 
and imaginary parts of the inverse Green's function at 
the same frequency) is satisfied nearly at 2J or 0.8t for 
the undoped case (Fig. ||(a)), just a few percent doping 
changes both the real and imaginary parts so drastically 
near 2 J that the quasiparticle condition is never met as 
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FIG. 2. Inverse of the interacting Green's functions at (0, 0) 
for two different doping concentrations. The solid and dotted 
curves denote the real and imaginary parts, respectively, (a) 
The doping concentration 0.0%. (b) The doping concentra- 
tion 5.54%. 



shown in Fig. @(b). On the other hand, the change of 
the self-energy at (7r/2,7r/2) (which is not shown here) 
is much smaller than that at the origin and thus the 
quasiparticle state persists at (tt/2, 7t/2). The disappear- 
ance of the quasiparticle state by doping is found in the 
significant portion of the Brillouin zone near the origin. 
This is mainly responsible for the dynamical narrowing 
of the quasiparticle band by doping, since the quasipar- 
ticle states near the origin comprise the large portion 
of the upper part of the quasiparticle band at x = 0.0. 
The total quasiparticle bandwidth is approximately 2J 
for 0.0% doping, while it is reduced to less than 0.8 J for 
5.54% doping by dynamical effects. We also found that 
doping not only makes high lying states near (0, 0) point 
unstable, but also strongly modifies low lying states. As 
the doping concentration is increased, the Fermi surfaces 
centered at (±7r/2, ±7r/2) with an ellipsoidal shape be- 
come so rapidly elongated along the antiferromagnetic 
zone boundary that they become a large (connected) 
Fermi surface about 10% doping level (Details will be 
published elsewhere ||). 

Figures ||(a)-(d) show Bn(q,ii) at the smallest mo- 
mentum (0, 7r/10) in the present calculation for several 
doping concentrations. In the undoped case spin wave ex- 
citations are eigenstates of our model Hamiltonian where 
they have an infinitely long lifetime as seen in Fig. [|(a). 
For 0.95, 2.93, 5.54 % doping concentrations, the renor- 
malization of the peak strength is small and we find only 
their small downward shifts, namely, softening of the spin 
excitation velocity. For 2.93 % doping (Fig. ||(c)), how- 
ever, a pronounced peak occurs near the origin, although 
its spectral weight is relatively smaller than that of the 
main peak. The occurrence of this peak can be under- 



FIG. 3. Spectral functions of the diagonal component of 
the matrix spin excitation propagator, Bn(q,v), at (0, 7r/10) 
and at (n/2, tt/2) for four different doping concentrations, (a), 
(b), (c) and (d) are Bu(q,v) at (0,tt/10) for 0.0, 0.95, 2.93 
and 5.54%, respectively, (e), (f), (g) and (h) are B\i(q,u) at 
(7r/2,7r/2) for 0.0, 0.95, 2.93 and 5.54%, respectively. 



stood more clearly by examining the real and imaginary 
parts of the inverse Dn(q,v), as shown in Figs. [I| For 
0.95 % concentration (Fig. |](a)), the pole condition is 
fulfilled only at 0.1 Qt and some small structure due to 
doping develops near the origin. For 2.93 % doping 
(Fig. ||(b)), this new structure grows so drastically that 
the pole condition is also satisfied near the origin. Its 
small spectral weight is caused by its small residue at that 
frequency. At 5.54 % doping (Fig. 0(c)), the real part 
grows further and the pole condition is less fully satisfied 
near the origin. Thus, the new structure near the origin 
becomes smaller for higher than 2.93 % doping concen- 
trations. Evidently the appearance of the new mode near 
the origin happens at between 0.95% and 2.93% doping 
within our numerical resolution. This corresponds to the 
disappearance of antiferromagnetic long-range order, as 
first pointed out by Sherman et al. |(| These authors ar- 
gued that the appearance of zero or negative frequency 
boson modes at zero temperature yields an infinite num- 
ber of bosonic excitations, implying the inadequacy of 
the spin wave basis approximation in that momentum. It 
indicates the disappearance of long-range antiferromag- 
netic order. This can be also gauged by the behavior of 
the anomalous (off-diagonal) spin excitation propagators. 
We found that the off-diagonal spin propagators become 
significant near the origin, (0,0), while they are vanish- 
ingly small near the zone boundary in spite of their strong 
renormalization and damping. This means the conden- 
sation of long wavelength spin excitation modes, just as 
electrons condensate in the superconducting state. Ap- 
parently this is an unphysical result, caused by the spin 
wave basis approximation for that particular momentum. 
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FIG. 4. Inverse of the diagonal component of the matrix 
spin excitation propagator at (0, n/10) for three different dop- 
ing concentrations. The solid and dotted curves denote the 
real and imaginary parts, respectively, (a) The doping con- 
centration 0.95%. (b) The doping concentration 2.93%. (c) 
The doping concentration 5.54%. 



The main spin excitation peak in Fig. ||(b) is slightly 
shifted downward for 0.95 % doping, leading to the soft- 
ening of the mode. We evaluated the concentration where 
the peak occurs at zero frequency by employing a sim- 
ple interpolation scheme to higher doping cases, which is 
essentially the same method as Igarashi et al. || used. 
These authors calculated the amount of softening (shift) 
of the mode by explicitly taking into account coherent 
and incoherent contributions and also by using pertur- 
bation theory to obtain the spectral function for a finite 
density of holes from that for a single hole. According 
to our simple interpolation, the complete softening of 
the main peak happens approximately at 11% in good 
agreement with 12% found by the above authors. This 
is evidently too high to be reconciled with experiments. 
The present calculation shows that the disappearance of 
long-range antiferromagnetic order is not determined by 
the complete softening of the mode but instead by the 
condition that the velocity of the long wavelength spin 
excitation becomes smaller than the Fermi velocity. As 
doping concentration is increased, the spin excitation ve- 
locity decreases and at the same time the Fermi veloc- 
ity rapidly increases as \[x in the small doping limit. 
Hence, self-consistent treatment of hole and spin dynam- 
ics is very important to predict a small value of the crit- 
ical concentration. The softening and damping of the 
spin excitation peaks at high momenta near the antifer- 
romagnetic zone boundary is quite drastic even for small 
doping, as exhibited in Figs, ^(e)-(h). In spite of this 
significant change in the spectral properties of spin exci- 
tations at these momenta, the off-diagonal (anomalous) 



components of the spin excitation propagator is vanish- 
ingly small and thus spin wave basis approximation ap- 
pears still valid, in agreement with the report by Sokol 
et al. 

In summary, the evolution of hole and spin dynamics 
in high T c superconductors was studied within the self- 
consistent noncrossing approximation in the small den- 
sity limit of doped holes. As the doping concentration is 
increased, string excitations rapidly disappear implying 
significant loss of long-range electron correlations. High 
energy states also lose their quasiparticle properties com- 
pletely and thus the quasiparticle energy band becomes 
considerably narrower than that of the undoped case. 
The strong change of the low lying quasiparticle band 
caused by dynamical effect of doped holes by exchange 
of spin excitations shows a tendency that separated small 
hole pockets become a large Fermi surface at relatively 
small doping concentrations. The disappearance of long- 
range antiferromagnetic order is found at between 0.95% 
and 2.93% doping levels both in the diagonal and off- 
diagonal components of the spin excitation propagator. 
The spin excitations at high momenta become strongly 
renormalized and damped by doping, yet they are still 
well described within spin wave basis approximation in 
agreement with the previous report. 
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